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rhe equation of state of cold baryonic matter is studied within a relativistic mean-field model with hadron masses and coupling 
constants depending on the scalar field. All hadron masses undergo a universal scaling, whereas the couplings are scaled differently, 
rhe appearance of hyperons in dense neutron star interiors is accounted for, however the equation of state remains sufhciently stiff 
f the reduction of the (f) meson mass is included. Our equation of state matches well the constraints known from analyses of the 
istrophysical data and particle production in heavy-ion collisions. 


1. Introduction 


A nuclear equation of state (EoS) is one of the key ingre¬ 
dients in the description of neutron star (NS) properties [1], 
mpernova explosions [2] and heavy-ion collisions [3,4]. A 
comparison of various EoSs in how well they satisfy var¬ 
ious empirical constraints was undertaken in Ref. [5] for 
the EoSs obtained within relativistic mean-field models 
(RMF) and some more microscopic calculations and in 
Ref. [6] for the Skyrme models. It turns out difficult to rec- 
Dncile the constraint on the maximum NS mass, which must 
be larger than 1.97 Mq after the recent measurements re¬ 
ported in [7,8], and the upper constraints on the stiffness of 
the EoS extracted from the analyses of heavy-ion collisions 
[HICs) [3,4]. Another relevant constraint on the EoS of the 
NS matter is imposed by the direct Urea (DU) processes, 
like n p -\- e v, which occur as soon as the nucleon 
density exceeds some critical value tiqu- The occurrence of 
these very efficient processes, even with account for the nu¬ 
cleon pairing, is hardly compatible with NS cooling data, 
if the value of the NS mass, at which the central density 
becomes larger than n^u, is < I-SMq (the so-called 

“strong” DU constraint) [9,5]. There should be Mgy < 
1.35M0 (the “weak” DU constraint) [10,5], since 1.35 Mq 
is the mean value of the NS mass distribution, as it follows 
from the analysis of the observational data on NSs in bi¬ 
nary systems. The DU problem appears in the EoSs with 
linear dependence of the symmetry energy except, maybe, 
most stiff ones. All the standard RMF EoSs and the mi¬ 
croscopic Dirac-Brueckner-Hartree-Fock (DBHF) EoS suf¬ 
fer of this linear dependence. On the contrary, variational 
calculations of the Urbana-Argonne group with A18 -I- 5v 


UIX* forces [11], as well as the RMF models with den¬ 
sity dependent hadron coupling constants [12], generate a 
weaker growth of the symmetry energy with the density, 
and the problem with the DU reactions is avoided. The later 
models are also able to describe NSs, as heavy as those in 
Refs. [7,8]. 

The problems worsen if strangeness is taken into account, 
because the population of new Fermi seas of hyperons leads 
to a softening of the EoS and reduction of the maximum 
NS mass. By employing a recently constructed hyperon- 
nucleon potential, the maximum masses of NSs with hy¬ 
perons are computed to be well below 1.4Mq [13]. Also, 
within RMF models one is able to explain observed massive 
NSs only, if one artificially prevents the appearance of hy¬ 
perons, cf. [13,14] and the references therein. This is called 
in the literature, the “hyperon puzzle”. So, the difference 
between NS masses with and without hyperons proves to 
be so large for reasonable hyperon fractions in the stan¬ 
dard RMF approach that in order to solve the puzzle one 
has to start with very stiff purely nuclear EoS, that hardly 
agrees with the results of the microscopically-based vari¬ 
ational EoS [11] and the EoS calculated with the help of 
the auxiliary field diffusion Monte Carlo method [15]. Such 
an EoS would also be incompatible with the restrictions 
on the EoS stiffness extracted from the analysis of nucleon 
and kaon flows in heavy-ion collisions [3,4]. All suggested 
explanations require additional assumptions, see discussion 
in [16]. For example, the inclusion of an interaction with a 
(/>-meson mean field, and the usage of smaller ratiohypreon- 
nucleon coupling constants following the SU(3) symme¬ 
try relations [17], as well as other modifications performed 
within the standard RMF approach, all help to increase the 




NS mass. 

There is also another part of the “hyperon puzzle”, which 
attracted less attention so far. With the hyperon coupling 
constants introduced with the help of the SU(6) symmetry 
relations the critical densities for the appearance of first 
hyperons prove to be rather low, n^u ~ 3no, cf. [18,19]. 
However with the appearance of the hyperons the efficient 
DU reaction on hyperons, e.g. A —>■ p + e + S', occurs that 
may potentially cause a very rapid cooling of the NSs with 
M > M]^u, being the NS mass, at which the central 
density reaches the value The problem should be addi¬ 
tionally studied with account for a suppression of hyperon 
concentrations and weak interaction vertices of hyperons 
compared to the corresponding neutron concentrations and 
nucleon weak interaction vertices, and with account for the 
hyperon pairing, cf. [20,21]. 

In Ref. [10] two of us formulated an RMF model, in 
which hadron masses and meson-baryon coupling con¬ 
stants are dependent on the a mean field. A working model 
MW(n.u., z=0.65) labeled in [5] as KVOR model has been 
constructed. This model was shown in [5] to satisfy appro¬ 
priately the majority of experimental constraints known 
by that time. In Ref. [22] the particle thermal excitations 
were incorporated and the model was successfully applied 
to description of heavy-ion collisions. However hyperons 
are not included in the model. Even without hyperons the 
KVOR EoS with the added (BPS) crust EoS from [23] 
yields = 2 .OIM 0 that fits the new constraint [7,8] 

only marginally. 

In the present Letter we will show that within the RMF 
models with hadron masses and coupling constants depen¬ 
dent on the a mean field one is able to overcome the men¬ 
tioned above problems and to construct the appropriate 
EoS with hyperons satisfying presently known experimen¬ 
tal constraints. 


2. Energy density functional 

In our RMF model we include nucleons N = (p, n) and 
hyperons H = (A°, interacting with mean fields 

of mesons m = a,u!, For simplicity, we drop the p me¬ 
son self-interaction and disregard, therefore, a possibility 
of the charged p-meson condensation discussed in [10]. As 
the bare masses of the particles we take ttin = 938 MeV, 
iTiA = 1116 MeV, ms = 1193 MeV, ms = 1318 MeV and 
neglect the small mass splitting in isospin multiplets. Lep¬ 
ton masses are rrie = 0.5 MeV and = 105 MeV. 

Following the approach [10], the scalar field enters as a di¬ 
mensionless variable / = ga-NX(TN{o’)a/mN and the meson- 
baryon coupling constants gmB, B = {N,H), are made 
(j-dependent with the help of scaling coupling constants, 

gmBXmBif) with Xu>H{f) = Xu^nU) , Xpuif) = XpnU) ■ 
The bare masses of hadrons ruB and rrim are replaced in 
the model by the effective masses = niB^Bif) and 
= rrim^mif) scaled by the functions 


^H{f) = ^ - XaH—^i<yHf , (1) 

mn 

where x^h = g^n/gaN and = XcrH{f)/XaN{f)- 

Taking into account the equations of motion for vector 
fields, the energy density of the cold infinite matter with 
an arbitrary isospin composition is recovered from the La- 
grangian of the model in the standard way, see Ref. [10]: 
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where X/^^p'^B guj{p)B/guj{p)N^ X(pB gpH/gujNt g^N 0. 

The values of the isospin projections for baryons follow 
from the Gell-Mann-Nishijima relation tsB = Qs — (1 + 
Sb)/ 2, where Qb and Sb are the baryon electric charge and 
strangeness, respectively. The baryon densities are related 
to the baryon Fermi momentum as ns = Pp and the 

fermionic kinetic energy density is 


Ekin(PF,m) 



p^dp 


XjP + rnP 


The dimensionless coupling constants are Cm = 
for all m except (j>. Since the ratios X(f,H are determined 
through gu:N, the 0-field contribution enters the energy 
density with the constant Ccf, = Cujmuj/mcj,. Here we take 
niu) = 783 MeV, = 1020 MeV. Bare masses and cou¬ 
pling constants of all mesons except </> enter the energy den¬ 
sity only in combinations Cm and the scaling functions d’m 
and Xm enter only through the scaling factors 


Vmif) = ‘^mU)/XmNif) ■ ( 3 ) 


Therefore we actually do not need to determine ^m{f) and 
Xm(f) separately, but only Pmif) combinations. The self¬ 
interaction of the scalar field introduced usually in RMF 
models through a potential U (/) is hidden now in the def¬ 
inition of Paif)- The equation of motion for the remaining 
field variable / follows from the minimization of the energy 
density dE[f]/df = 0 . If we suppress </> and H terms, put 

Per = 1 -k 2^(|/^ -k |/^) and put all other scaling func¬ 
tions to unity, we recover the energy density functional of 
the standard non-linear Walecka a-to-p model. 

For the NS matter sustained in the /3-equilibrium the 
Fermi momenta of a baryon can be expressed through the 
baryon chemical potential, PB) as pp g = {pb—Vb)^ — 

m|$|, where UBwir = Clx^BriB/Pi^ + ClxpBhBni/pp + 

C^XfftBns/P4>] Pb is related to the nucleon and electron 
chemical potentials as pB = Pn — QbPb- Solving the sys¬ 
tem of equations for pf,b and making use of the electro¬ 
neutrality condition ^b QbBb = ne + n^, where the lep¬ 
ton densities are Ug = = {pi — m^)^/^/37r^, 

we can express the hadron densities and the total energy 



density through the total baryon density n = ub- The 
total pressure is calculated as P = — E. 

The parameters of the nucleon sector are tuned to re¬ 
produce the properties of nuclear matter at saturation: the 
saturation density no, the binding energy per nucleon Sq, 
the effective nucleon mass rn*jq, the compressibility modu¬ 
lus K and the symmetry energy J. The coupling constants 
of hyperons with vector mesons are interrelated by SU( 6 ) 
symmetry relations [26]: 

9ujK HcjS '^9u)N j 9‘^9^9pN 5 

2g0A = 2^02 = 9<t>=. = -^^9ujN , 9pA = 94,N = 0. (4) 

The coupling constants of hyperons with the scalar mean 
Held are constrained with the help of the hyperon binding 
energies per nucleon in the isospin symmetric matter 
(ISM) at n = no given by [18]: 

'^Wnd(’^o) = - {itin - rn*t^{no))xaH , (5) 

where we suppose ^a-Hifino)) = 1 and use 

^bind(^o) = -28 MeV, fbi„d(»^o) = 30 MeV, 
4lndK) = -15MeV. ( 6 ) 

The repulsive E potential prevents the appearance of E 
hyperons in all models considered below. 

The NS configuration follows from the solution of the 
Tolman-Oppenheimer-Volkoff equation. For n $ (0.6 — 
0.8)no, our RMF EoSs should be matched with the EoS 
of the NS crust, where the formation of a pasta phase is 
explicitly included. The presence of the pasta, although it 
may change transport properties of the matter, affects an 
EoS only slightly [24]. Therefore simplifying consideration 
we chose the frequently-used crust BPS EoS from [23], ig¬ 
noring the pasta phase. The same BPS crust EoS was used 
in [5] where it was joined with the various EoSs describ¬ 
ing the interior region. The pressure as a function of the 
density for the EoSs, we consider here, intersects the BPS 
pressure at density ni, such that 0.45no ^ ni ^ 0.7no- We 
construct the resulting pressure as a function of the den¬ 
sity by using the BPS pressure for n < 0.45no, then a cu¬ 
bic spline interpolation for 0.45no < n < O.lriQ, and the 
pressure for beta-equilibrium matter (BEM) given by our 
model for n > 0.7 ng. The choice of the interval, which in¬ 
cludes the intersection points for both EoSs, guarantees the 
smoothness of the interpolation. We have also checked that 
a narrowing of the interpolation interval limits almost does 
not reflect on such observables as the star mass and radius. 

3. Models 

We discuss now two models constructed according to the 
principles described above. One is the formal extension of 
the KVOR model from Ref. [10] to hyperons, which we call 
the KVORH model now. In this example we demonstrate 
the problems, which appear when one includes hyperons. 
Eor the second model, labeled as MKVOR (or MKVORH 
when hyperons are included), we propose new set of scaling 


Table 1 

Characteristics of KVOR and MKVOR models at saturation 


EoS ■ 

So 

no 

K 

m*^{no) 

J 

L 

K' 

-^sym 

[MeV] [fm-3] 

[MeV] 

[rriN] 

[MeV] [MeV] [MeV] [MeV] 

KVOR 

-16 

0.16 

275 

0.805 

32 

71 

422 

-86 

MKVOR 

-16 

0.16 

240 

0.73 

30 

41 

557 

-158 


functions. In Table 1 we present the saturation parameters 
for both models and the coefficients of the expansion of 
the nucleon binding energy per nucleon near the nuclear 
saturation density no, 

f = fo + y + P^j{n) + , 

J(n) = J + Le+^e^ + ... , (7) 

in terms of small e = {n — no)/3no and (3 = (n„ — np)ln 
parameters. 


3.1. KVORH model 

Reference [10] proposed a set of input parameters for 
the RMF model, which matches the APR EoS (in the rel¬ 
ativistic HHJ parameterization of [27]) up to n $ 4no, see 
Eq. (61) in [10]. To fulfill the DU constraint Ref. [10] intro¬ 
duced the scaling function r]p{f) ^ 1, see Eq. (63) in [10]. 
Thus the model labeled as KVR in Ref. [5] was constructed. 
The idea behind the KVOR modification of the KVR model 
was to demonstrate that introducing the additional scaling 
function 1 one can increase the maximum value of the 
NS mass without a sizable change of the EoS for densities 
n 4no. The coupling constants for the KVOR model are 
given in Eq. (58) in Ref. [10]. The KVOR model produces 
the maximum NS mass M^ax = 2 .OIM 0 , and the critical 
proton density for the DU reaction threshold on neutrons 
n^u = 3.92no corresponding to = 1.76M0. 

In the KVORH model the parameters x^h deduced from 
hyperon binding energies in Eq. (5) are 

= 0.599 , = 0.282 , = 0.305 . (8) 

The baryon concentrations for the KVORH model are 
depicted in Eig. 1 as functions of the baryon density n. 
For n < 0.6no the curves presented in Fig. 1 should be re¬ 
placed by those computed within the EoS of the crust. The 
KVORH model produces the maximum NS mass = 

1.66M0. The critical density for the appearance of first hy¬ 
perons, which is simultaneously the critical density for the 
onset of the DU reactions on hyperons (A in this case) is 
n^u = 2.82no, and the corresponding NS mass at which 
first As appear in the NS center is = 1.38M0. The 
total strangeness concentration (the ratio of the number of 
strange quarks to the total number of quarks) in the NS 
with the maximum mass is fs = 0.034. 




n/n„ n/n„ 

0 0 


Fig. 1. Baryon concentrations and the variable / in BEM as func¬ 
tions of total baryon density for the KVORH, MKVORH^ and 
MKVORH^cr models. Dotted lines show the thresholds of the DU 
reaction on neutrons. For the MKVORH0 and MKVORH^cr models 
the thresholds are not distinguishable on the plot scale. Short ver¬ 
tical bars on the lines show the maximum densities reachable in the 
NS for the given model. 


Table 2 

Parameters of the MKVOR model 


Cl 

Cl 
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c- 10^ 

d 

a 

z 

Cl(jj 

234.15 

134.88 81.842 

4.6750 

-2.9742 

-0.5 

0.4 

0.65 0.11 

fuj 

^ 1 fp 

Clp 


4=^ 

4^' 

dp 

ep 

7.1 0.9 

3.11 28.4 0.522 

0.448 

-0.614 

3 

0.8 

-4 

6 


in the DD and DD-F models with the density dependent 
coupling constants and bare meson field masses, cf. [25,5]. 

The maximum NS mass increases in the MKVOR model 
compared to the KVOR model up to Mmax = 2.33 M©. 
The DU threshold values are ripy = 4.14no and Mgy = 
2.22Mq. 

In the model with hyperons (MKVORH) we use the aH 
coupling constant ratios deduced from hyperon binding en¬ 
ergies in ISM following Eq. (5) for = 1; 

a;^A= 0.607, = 0.378 , 0.307. (10) 


3.2. MKVORH model 


Reference [10] showed that the EoS is more sensitive to 
the value of m*^(no) than to the compressibility K. The 
smaller m^(no) is in a certain RMF model, the larger is the 
value of the maximum NS mass. The input parameters for 
the new MKVOR model are listed in Table 1 together with 
the corresponding parameters of the nuclear binding energy 
per nucleon at saturation. We took in the MKVOR model 
a smaller value of m^(no) than in the KVOR model and 
a smaller value of the compressibility, K = 240 MeV, that 
agrees with canonical value K = 240 ± 20 MeV extracted 
from the analysis of giant monopole resonances (GMR) [28]. 

The scaling functions of the MKVOR model are as fol¬ 
lows: 


^U(/) = 1 - lclbf[i -f J(c - lclb^)f] + Idf , 

Vu>if) = ( X+1/° ) + + ’ 

Vpif) = + ''' ({) 

1 + Op 


+13 exp 


(/-/p)^(l + ep(/-/o)^) ] 

^l + dp(/-/o) + ep(/-/o)2J ■ 


(9) 


The parameters of the model and of the scaling functions 
are collected in Table 2. The choice of the scaling functions 
is definitively not unique. Their final form was tuned to 
satisfy best the experimental constraints, that we demon¬ 
strate below, and to keep a connection to the KVOR model 
parameterization. The first term in is the same as in the 
KVOR model, the function rja- and the first three terms in 
rjp are basically the reparameterization of the functions of 
the KVOR model. The new terms, the second one in 
and the last one in rjp are added to control the growth of 
the scalar field with density. Note that the density depen¬ 
dencies of our scaling functions for the aN, ojN and pN 
coupling constants prove to be similar to those exploited 


Note that the value x^a — 0.61 is close to the best value 
derived from hypernuclei XaA — 0.62 in Ref. [26]. 

To have an opportunity for an increase of the maximum 
NS mass in the model with hyperons, we incorporate the 
0 -meson mean field with a scaled 0 meson mass (version 
labeled MKVORH0) and, additionally, allow for a scal¬ 
ing of the aH coupling constants, ^aH(n) ^ 1 (version 
MKVORH0cr). 

In the model MKVORH0 we will exploit the very same 
scaling of the 0 -meson mass as for other hadrons $0 = 
I — /. This implies the scaling function 

W = (l-/)^- (11) 

We use the minimal model, assuming parameterization 
Eq. (4) for the vector-meson-hyperon coupling constants, 
and the aH coupling constants from Eq. (10). As the re¬ 
sult the maximum NS mass in MKVORH0 model becomes 
Afmax = 2.22 Mq with the strangeness concentration 
/s = 0.023. The critical density for the appearance of first 
hyperons is = 2.63no, corresponding to the star mass 
Mpu = 1.43 M©. Although the model does not satisfy 
the “strong” DU constraint (Mdu > 1-5M©), we should 
notice that for reactions on hyperons this constraint might 
be a bit soften since the baryon part of the squared ma¬ 
trix element of the DU reaction on As is 25 times smaller 
than that for the DU reaction on neutrons [20], besides the 
pairing gaps might be not as small. 

The effect of the (n) scaling we demonstrate at hand 
of the MKVORH0 (t model, where we take (n) such that 

< "Uo) = 1 , and assume that ianiji) decreases with 
an increase of n and vanishes for densities n > minlnpy}. 
This means that effectively we will exploit vacuum hyperon 
masses for n > min{n((u}. Note that the KVOR model ex¬ 
tended to high temperatures in Ref. [22] (called there as the 
SHMC model) matches well the lattice data up to T ^ 250 
MeV provided aB coupling constants for all baryons ex¬ 
cept the nucleons, are artificially suppressed, that partially 
motivates our choice of suppressed values 

























Fig. 2. (Left) The symmetry energy coefficient J(n) as a function of 
the nucleon density calculated in the KVOR and MKVOR models. 
Shaded area shows the constraint from the study of analog isobar 
states (AIS) in [29]. Hatched area is the constraint from the electric 
dipole polarizability (ao) in ^^®Pb [30]. (Right) The nucleon optical 
potential as a function of the nucleon kinetic energy for ISM at n = no 
for the KVOR, cf. [22], and MKVOR models. Shaded area shows the 
extrapolation from finite nuclei to the nuclear matter from [32]. 

The baryon concentrations from the MKVORHc/) and 
MKVORH^cr models are shown in Fig. 1. The proton 
fractions of the MKVORHc^ and MKVORH(/)cr models are 
smaller than those for the KVORH model. We also see 
that inclusion of the (j) scaling (11) reduced the hyperon 
population. The reduction of the aH coupling constants 
prevents the appearance of A and hyperons and shifts 
the threshold density of the appearance to higher val¬ 
ues. Without As the reaction 2“ —>• K + e~ + Df. does 
not occur and the DU threshold is determined by the DU 
reactions on nucleons. Replacing the values of f(n) de¬ 
picted in Fig. 1 for the BEM in Eq. (1) one can recover 
the density dependence of the effective hadron masses and 
from Eqs. (9) and (11) that of the effective coupling con¬ 
stants. Within the MKVORH(()cr model we get Mmax = 
2.29 Mq, n^u = 3.69no (Mgy = 2.09 Mq), and the total 
strangeness concentration in the heaviest NS is reduced to 
/s = 6.2 X 10-3. 

Applying the (/)-mass scaling and the scaling to the 
KVORH model we obtain for the KVORH^ model M^ax = 
1.88 Mq and that the first hyperons, As, appear at the den¬ 
sity = 2.81 no (M^^ = 1.37 Mq). The strangeness 
fraction is /s = 0.035. For the KVORH(()cr model we find 
Almax = 1-96 Mq, fs = 9.2 X 10-3. The first among hyper¬ 
ons appear 2-s, therefore the DU threshold is shifted to 
nBu = 3.95 no (Mgy = 1.76 Mq). 

Below we compare how well the EoSs obtained within 
the MKVORH^cr and KVORH models satisfy various phe¬ 
nomenological constraints. 

4. Constraints on the models 

4.1. Symmetry energy and nucleon optical potential 

Constraints on the density dependence of the symmetry 
energy [parameter J(n) in Eq. (7)] are extracted in [29] 
from the study of the analog isobar states and in [30] from 


the electric dipole polarizability of ^°®Pb nuclei. They are 
shown in Fig. 2 (left panel) by the shaded and hatched re¬ 
gions, respectively, together with the symmetry energies 
calculated in the KVOR and MKVOR models. We see that 
the both models follow the lower boundary of the region. 

The dependence of the nucleon optical potential on the 
nucleon kinetic energy in the ISM at n = no is shown in 
Fig. 2 (right panel). The shaded region is extracted from 
the atomic nucleus data [31] and recalculated to the case of 
the infinite nuclear ISM in [32]. The KVOR model describes 
the nucleon optical potential for low and high particle en¬ 
ergies but does not describe it for intermediate energies. 
The MKVOR model describes the nucleon optical poten¬ 
tial rather well for nucleon energies Em — ruM ~ 400 MeV. 
To fit appropriately the data at higher particle energies, 
the momentum dependence of the NN interaction would 
be required that is not present in the mean-field approach. 
The iso-vector part of the optical potential is 

less constrained by the data, therefore we do not show it. 

4.2. Particle flow in heavy-ion collisions 

The analysis of the transverse and elliptical flow data in 
HICs allowed [3] to extract a constraint on the pressure of 
the ISM as a function of the nucleon density and to recon¬ 
struct the pressure for the purely neutron matter (PNM) 
with some assumptions about the density dependence of 
the symmetry energy J{n) (soft or stiff one). Reference [33] 
on the basis of calculations of the kaon production in HICs 
[4] provided some restrictions on the pressure at lower den¬ 
sities. These constraints are shown in Fig. 3 together with 
the results of the GMR data analysis taken from [33] and 
the pressure calculated in the KVOR and MKVOR mod¬ 
els. The constraints rule out a very stiff EoS. We see that 
the KVOR model satisfies the requirements. The MKVOR 
model fulfills the constraints, for n < 4no in ISM. For PNM 
the MKVOR curve passes through the hatched region with 
stiffer J(n), whereas the KVOR EoS is softer. The choice of 
the smaller TO^(no) leads to a stiffening of the EoS in ISM 
and the scaling functions rja-if) and rjuiif) are chosen to 
soften the EoS for n 4no to fulfill the nucleon- and kaon- 
flow constraints. To increase the maximum NS mass we had 
to stiffen the EoS in the BEM that was accomplished by 
the choice of the rjp function. 

4.3. DU constraint 

In the BEM the DU process on neutrons, n p-\-e-\-V'e, 
can occur only, if the proton fraction is high enough so 
that the Fermi momenta of neutrons, protons and electrons 
iPF,^=n,p,e) satisfy the inequality pp.n < PF,p + PF,e- Usu¬ 
ally, RMF models yield uncomfortably low values of the 
threshold densities for these reactions and correspondingly 
low values of the NS mass, M^^, at which the process be¬ 
gins to occur in the NS center. Every star with a mass only 
slightly above Mgy cools down fast due to the DU process. 





Fig. 3. Pressure as a function of the nucleon density in ISM and 
PNM for the KVOR and MKVOR models. Hatched areas show the 
empirical constraints from the analyses of a particle flow in HICs 
in [3], kaon production in HICs [33,4], and the GMR data [33]. 

even in the presence of nucleon pairing, and becomes invis¬ 
ible for the thermal detection within few years [9]. Most of 
single NSs have likely masses below I.SMq in accordance 
with the type-II supernova explosion scenario [2] and the 
population synthesis analysis [34]. Therefore, it is natural 
to believe that majority of the pulsars, which surface tem¬ 
peratures have been measured, have masses M < 1.5Mq. 
The analysis of these data in the existing cooling scenarios 
supports the constraint Mgy > 1.5Mq [35,36]. The ade¬ 
quate description of the new data on the cooling of the Cas- 
siopea A also requires the absence of the DU reactions [35]. 

In the presence of A hyperons, reactions K ^ p + e -\- D 
may occur. As seen in Fig. 1, proton concentrations for the 
KVORH, MKVORH^ models do not exceed the neutron- 
DU threshold for n < riQUi ^or higher densities the 
DU reactions on As start occurring. For the MKVORHi^cr 
model As do not appear at all and the DU processes occur 
for n > npy. 

4.4. Gravitational mass versus haryon mass constraint 

The unique double-neutron-star system J0737-3039 with 
two millisecond pulsars provided an important constraint 
on the nuclear EoS. The gravitational mass of one of the 
companions (B) is very low Mq = 1.249 ± O.OOIMq [37] 
which implies a very peculiar mechanism of its creation 
- a type-I supernova of an 0-Ne-Mg white dwarf driven 
hydrostatically unstable by electron captures onto Mg 
and Ne. Knowing this mechanism Refs. [38,39] calculated 
the number of baryons in the pulsar and the correspond¬ 
ing baryon mass: Mb = 1.366-1.375M0 [38] and Mb = 
1.358-1.362 [39]. The constraint of Ref. [38] can be released 
by 1% Mq because of a possible baryon loss and a critical 
mass variation due to carbon flashes during the collapse. 
Therefore one can speak about “strong” (without the mass 
loss) and “weak” (with the mass loss) constraints on the 
EoS, respectively. Microscopically motivated EoSs, like the 
relativistic DBHF EoS [41], the APR EoS [11], the diffu¬ 
sion Monte-Carlo one [15], and many RME-based models 
do not fulfill the strong constraint of Ref. [38]. Many EoSs 



Fig. 4. (Left) Gravitational NS mass Mq versus the baryon mass 
Mq. Shaded rectangles show the constraints from J0737-3039(B) 
pulsar derived in Refs. [38] and [39] labeled by 1 and 2, respectively. 
The empty rectangle demonstrates the change of the constraint [38] 
by 1% Mq. (Right) The NS mass-radius relations for the KVORH, 
MKVORH^, MKVORH^cr models compared with constraints from 
the isolated NS RX J1856.53754 [40], QPOs in the LMXBs 4U 
0614-1-09 [42], the millisecond pulsar PSR J0437-4715 [43], and 
the Bayesian probability distribution analyses (BPA) [1]. The hor¬ 
izontal lines border the uncertainty range for the mass of PSR 
J0348-h0432 [8]. 

do not satisfy the constraint of [39] and even the weak 
constraint of Ref. [38], cf. [5]. 

In Eig. 4 (left panel) we plot the gravitational NS 
mass Me versus the baryon mass Mb- The KVOR model 
matches marginally the weak constraint of Ref. [38], 
whereas the MKVOR model matches marginally both the 
result from Ref. [39] and the strong constraint of [38], 
the latter was not reproduced by the EoSs considered in 
Ref. [5]. Note that hyperons do not appear in the NS with 
the mass of 1.25M0, therefore the curves for the models 
KVORH and KVOR coincide, as well as the curves for the 
MKVOR, MKVORH((i and MKVORH(/)(t models. Com¬ 
paring particle concentration for different models shown 
in Fig. 1 and the baryon mass of the NS, we observe a 
correlation: the smaller is the proton fraction within the 
density interval uq < n < 2.5no, the better the given EoS 
satisfies the baryon-mass constraint. Eor n ^ uq the value 
of the proton fraction is correlated with the values of J and 
L in Eq. (7). The value of J may vary only a little (from 
28MeV to 34MeV or even in a narrower interval) but L 
varies broadly in various works, e.g., cf. Fig. 4 in Ref. [33]. 
With a decrease of L curves in Fig. 4 (left panel) are shifted 
to the right. For L < 40MeV (at fixed other parameters) 
the MKVOR curves would pass through the shaded area 
1. The smaller L is for the given EoS, the less the proton 
fraction is in the relevant density interval no < n < 2.5no, 
and the better the constraint is satisfied. 

4.5. Mass and radius constraints 

In Fig. 4 (right panel) we show mass-radius relations of 
NSs for the KVORH, MKVORHcji and MKVORH(()cr mod¬ 
els. Largest precisely-measured masses of NSs are 1.97 ± 
O.O4M0 for PSR J1614-2230 [7] and 2.O1±O.O4M0 for PSR 
J0348-f0432 [8]. The MKVORH<() and MKVORHcjicr mod- 














els can describe these high-mass NSs, whereas the KVORH 
model fails badly. Experimental information about heavier 
NSs is plagued with large experimental errors and with ad¬ 
ditional theoretical uncertainties, see Ref. [1] for a review. 
In Fig. 4 (right panel) we confront our models with other 
constraints derived from the quasi-periodic oscillations in 
the low-mass X-ray binary 4U 0614-1-091 [42] and thermal 
spectra of the nearby isolated NS RX J1856.5-3754 [40]. 
More details on these constraints can be found in Ref. [5]. 
In contrast to the mass determination, there are no accu¬ 
rate estimates of NS radii. Some constraints were derived 
recently from the X-ray spectroscopy of PSR J04374715 
with the proper account for the system geometry [43] and 
from the Bayesian probability analysis of several X-ray 
burst sources in Ref. [1]. These constraints are also shown 
in Fig. 4. We see that the MKVORH(() and MKVORHc^cr 
satisfy the mass-radius constraints and produce the radii 
of NSs in a narrow interval 11.7 ± 0.5 km for star masses 
M > 0.5Mq. 

5. Concluding remarks 

We constructed relativistic mean-field models with 
scaled hadron masses and coupling constants including a 
(j) meson mean field and hyperons. The hyperon-vector- 
meson (w, p, </>) coupling constants obey the SU(6) symme¬ 
try relations. The most challenging is to fulfill the flow con¬ 
straint and produce a high maximum neutron-star mass si¬ 
multaneously. For that we introduced the scaling functions 
such that our equation of state is rather soft for n ^ 4no 
in the isospin symmetric matter but is sufficiently stiff in 
the beta-equilibrium matter. This behavior is achieved by 
the proper selection of the scaling functions tJuj and Pp. 
The inclusion of the (j) meson with the mass scaled in the 
same way as masses of nucleons and other mesons, but 
with the (j) coupling constants being fixed, allows to fulfill 
the empirical constraints on the maximum neutron star 
mass; see curve for the MKVORH(/) model in Fig. 4 (right 
panel). In such an approach the “hyperon puzzle” can be 
resolved within our model without any additional assump¬ 
tions, like, e.g., the change of a scheme for the choice of 
the hyperon-vector-meson coupling constants from SU(6) 
to SU(3), as in [17]. We stress that we would not succeed 
if we used the bare (j) meson mass. Other constraints are 
also satisfied except that the model MKVORH(() produces 
a rather low threshold value of the neutron star mass 
Mpu 1.44M0 for the occurrence of direct Urea reactions 
on As. On the other hand, such a value of might be 
already sufficiently high to cause no problems with the too 
rapid cooling of neutron stars, since the DU process on As 
might be less efficient than that on neutrons. Nevertheless 
we demonstrated how one can fully eliminate this possible 
deficiency. The model MKVORHt/xr, where the hyperon 
masses do not change in medium, satisfies appropriately 
all constraints discussed in this Letter, which are known 
from the analyses of atomic nuclei, heavy-ion collisions and 


neutron star data. Moreover, the maximum neutron star 
mass increased. As an interesting finding, we indicate that 
the smaller the proton fraction is in the density interval 
no < n < 2.5no and the smaller the value of L is at no, the 
better the baryon-gravitational mass constraint is fulfilled. 
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